Evolution of the Fermi surface of BaFe 2 (Asi_ ;! ;P x )2 
on entering the superconducting dome 
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Using the de Haas-van Alphen effect we have measured the evolution of the Fermi surface of 
BaFe2(Asi_ x P I )2 as function of isoelectric substitution (As/P) for 0.41 < x < 1 (T c up to 25 K). 
We find that the volume of electron and hole Fermi surfaces shrink linearly with decreasing x. This 
shrinking is accompanied by a strong increase in the quasiparticle effective mass as x is tuned toward 
the maximum T c . It is likely that these trends originate from the many-body interaction which gives 
rise to superconductivity, rather than the underlying one-electron bandstructure. 



Superconductivity in the 122 iron-pnictide family 
XFe2As2 (where X=Ba,Sr or Ca), can be induced by a 
variety of means, including doping [H, Q, pressure 0] or 
isoelectric substitution either on the Fe 4] or As sites 
5]. The highest T c achievable by each of these routes 
is roughly the same. It has been suggested [ESQ 
that the interband coupling between the hole and elec- 
tron sheets plays an important role in determining the 
magnetic or superconducting order formed at low tem- 
perature. Discovering how the Fermi surface evolves as a 
function of the various material parameters which drive 
the material from an antiferromagnetic spin density wave 
state, through the superconducting dome and eventu- 
ally towards a paramagnetic non-superconducting metal, 
should therefore be an important step toward gaining 
a complete understanding of the mechanism that drives 
high temperature superconductivity in these materials. 

In the antiferromagnetic state is it expected that the 
Fermi surface suffers a major reconstruction. This is 
supported experimentally by the observation [^ [Tot 11 1 



of very low frequency quantum oscillations in the un- 
doped XFe2As2 compounds, corresponding to very small 
Fermi surface pockets which are 1-2% of the total Bril- 
louin zone planar area. At the other extreme of the phase 
diagram, where the materials are paramagnetic and non- 
superconducting, quantum oscillation measurements of 
SrFe 2 P 2 and CaFe 2 P 2 show that the Fermi sur- 
face is in good agreement with conventional bandstruc- 
ture calculations. Up to now, tracking the changes in 
the Fermi surface across the phase diagram using quan- 
tum oscillations has not been possible because of the 
additional disorder and high introduced by doping. 
Measurements on the low T c (6 K) superconducting iron- 
pnictide LaFePO 14, lj| established that the Fermi sur- 
face is in broad agreement with bandstructure with mod- 



erate correlation enhancements of the effective mass. It is 
not clear whether the higher T c pnictide superconductors, 
which unlike LaFePO occur in close proximity to a mag- 
netic phase, are also well described by bandstructure and 
whether the electronic correlations change significantly. 

The substitution of P for As in the series 
BaFe2(Asi_ K Pj;)2 offers an elegant way to suppress mag- 
netism and induce superconductivity without doping 
As P and As are isoelectric, there is no net change 
in the ratio of electrons to holes and the system remains 
a compensated metal for different x (equal volumes for 
the electron and hole Fermi surfaces). This series has 
several remarkable properties which are similar to those 
observed in cuprate superconductors. Firstly, the tem- 
perature dependence of the resistivity changes from a 
quadratic (p ~ T 2 ) to linear behavior (p ~ T) as the 
system evolves from a conventional Fermi liquid (x = 1) 
towards the maximum T c (x = 0.33). Secondly, there 
is strong evidence from magnetic penetration depth [ijj , 
thermal conductivity 17 1 and NMR [llj measurements 
that for x = 0.33 the superconducting gap has line nodes. 

In this Letter, we report the observation of quantum 
oscillation signals in samples of BaFe2(Asi_ a: P 2 ;)2 as x 
is varied across the superconducting dome from x = 1 
to x = 0.41 with T c ~ 25 K (~ 0.8T c max ). Our data 
show that the Fermi surface shrinks and the quasiparti- 
cle masses become heavier as the material is tuned to- 
ward the magnetic order phase boundary, at which T c 
reaches its maximum. This implies that the significant 
change of electronic structure caused by many-body in- 
teractions plays a key role for the occurrence of uncon- 
ventional high-T c superconductivity. 

Single crystal samples of BaFe2(Asi- a; P 2 ;)2 were grown 
as described in Ref. [16j . The x values were determined 
by an energy dispersive X-ray analyzer. Measurements of 
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FIG. 1: (color online), a) Raw torque data for BaFe2(Asi_ a ;P :c )2 with x = 0.33, 0.41, 0.56 measured at 1.5 K. b) The oscillatory 
part of torque for selected samples and c) the corresponding Fourier transform for x — 0.41, 0.56 between 35-55 T at 1.5 K and 
x = 0.64, 0.71 between 20-30 T at 0.35 K as well as x = 1 between 12-16.8 T at 0.092 K. The inset shows the oscillatory signal 
for x = 1.0. 



the quantum oscillations in magnetic torque [the de Haas- 
van Alphcn (dHvA) effect] were made using a miniature 
piezo-resistive cantilever technique. Experiments were 
performed in: a dilution refrigerator system with DC 
fields up to 17 T (Osaka), a pumped 3 He system with 
DC fields up to 30 T (Nijmegen) and 45 T (Tallahas- 
see) and a pumped 4 He system with pulsed fields up to 
55 T (Toulouse). Our data are compared to band struc- 
ture calculations which were performed using the wien2k 
package jl9j |. 

Fig- da) shows raw magnetotorque t(H) data for three 
values of x measured up to 55 T at T ~ 1.5 K, with the 
magnetic field direction close to the c-axis. The torque 
response for the higher T c samples (x = 0.33, T c ~ 30 K 
and x = 0.41, T c ~ 25 K) is highly hysteretic. [ioHi rr 
increased substantially with x reaching a maximum of 
51.5 T for the highest T c sample (x = 0.33), which is 
close to the estimated H c i value [17|. By subtracting 



a smooth polynomial background the oscillatory dHvA 
signal is clearly seen above H irr (Fig.QJb)) for all samples 
except for that with the highest T c (x = 0.33). 

From the fast Fourier transform (FFT) spectra of the 
oscillatory data (Fig.QJc)) we can extract the dHvA fre- 
quencies F. These are related to the extremal cross- 
sectional areas, Ak, of the Fermi surface orbits giving 
rise to the oscillations via the Onsager relation, F — 
(TL/2ire)Ak- From the evolution of these frequencies (see 
Fig. [2]) as the magnetic field is rotated from being along 
the c axis (9 = 0°) towards being perpendicular the c 
axis (9 — 90°), we can deduce the shape of the Fermi 



surface. 

Band structure calculations of the Fermi surface of the 
end members of the BaFe2(Asi_ :r P :E )2 series are shown in 
Fig. [2] For both BaFe2 As2 and BaFe2P2 the calculations 
were done with the experimental lattice parameters and 
pnictide z position in the non-magnetic state. The two 
electron sheets at the zone corner are quite similar in size 
and shape in both compounds but there are significant 
differences between the hole sheets. In BaFe2As2 three 
concentric quasi-two-dimensional hole tubes are located 
at the zone center, whereas in BaFe2P2 the inner one of 
these tubes is absent whereas the outer tube has become 
extremely warped. The Fermi surface of BaFe2P2 is quite 
similar to that for SrFe2P2 (l2| . 

For SrFe 2 P 2 (i2|, CaFe 2 P 2 and LaFePO it was 
observed that the strongest dHvA signals came from the 
electron sheets and so this is likely also to be the case 
for BaFe2P2- Thus, we assign the two strongest peaks (a 
and 0) to the inner and outer electron sheet, respectively 
(see Fig. [He)). The observed frequencies for x = 1 are 
somewhat smaller than those predicted by the calcula- 
tions and a rigid band shift of 45 — 75 meV is needed to 
bring the orbits into agreement with experiment. This 
is quite comparable to the shifts 50 — 60 meV needed for 
the electron sheets in SrFe 2 P2 [13] • Although there are 
other weak peaks in the FFT which could come from the 
hole sheets, further measurements with better signal to 
noise ratio are needed to confirm these. 

As the sample composition is varied towards optimal 
doping, the signal from the electron sheets is reduced but 
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FIG. 2: (color online), a) Angle dependence of the predicted orbits from B ||[001] (0 = 0°) to B ||[1OO](0 = 90°) for x = 1 
(BaFe2P2) and the electron orbits observed experimentally for b) x = 1, c) x = 0.72, d) x = 0.64, e) x = 0.56. Solid lines 
correspond to F(0 = 0)/cos#. The insets show the temperature dependence of the Fourier amplitudes for each composition, 
which are fitted (solid line) to the Lifshitz-Kosevich formula Xj sinhX, X = 14.69 m*T/B to determine the effective masses, 
m* [2C| . The calculated Fermi surfaces of the end members, f) BaFe2As2 (non-magnetic) and g) BaFe2P2 are also shown. 



two peaks are clearly visible over the full doping range in 
sufficiently high fields (see Fig.QJc)); the mean free path 
for the (3 orbits decreases from I ~ 800 A to 170 A when 
x varies from x = 1 to 0.41. Importantly, we observe 
that the frequency of both these electron orbits decreases 
linearly with decreasing x (see Fig. EJa)). For x = 1 the 
(3 frequency equals 2.30 kT which decreases to 1.55 kT 
for x — 0.41. Although we do not clearly observe them, 
charge neutrality means that the hole sheets must shrink 
by the same factor. 

The effective masses of the quasiparticles on the var- 
ious orbits were determined by fitting the temperature 
dependent amplitude of the oscillations to the Lifshitz- 
Kosevich formula [20| (see insets of Fig. [2|). As shown 
in Fig. EJb), the effective mass increases significantly as 
x approaches the spin density wave ordered phase. At 
the same time T c increases and shows a maximum at the 
boundary of the magnetic order. As no such increase is 
expected from the bandstructure (see later) the result im- 
plies a significant rise in the strength of the many-body 
interactions, most likely caused by spin-fluctuations. A 
similar striking rise in the mass enhancement close to an 
antiferromagnetic quantum critical point has also been 
observed in quasi-two dimensional heavy-Fermion sys- 
tems [Hj]. 

For x = 1 bare masses calculated from the bandstruc- 
ture for the electron sheets vary from 0.8 to 1.0 m e . 
To obtain a dHvA frequency of 1.55 kT, appropriate for 
x = 0.41, requires a rigid shift of the energy of band 4 
(see Fig.[2l[a)) by 170 meV. The bare mass then decreases 
from 0.92 to 0.81 m e , in sharp contrast to the large rise 
seen in experiment. The many body enhancement factor 



m* /mi) therefore increases from ~ 2 to ~ 4 as x decreases 
from 1 to 0.41. Note that conventional electron-phonon 
coupling is weak in these materials, and is estimated to 
only account for ^25% of the observed mass enhance- 
ment [II]. 

The shrinking of the electron and hole sheets as x 
decreases could come from either the underlying one- 
electron bandstructure or alternatively it could be driven 
by many-body correlation effects. We have estimated the 
change in the electron sheet area by calculating the band- 
structure for both BaFe2As2 and BaFe2P2 with values of 
the lattice constants and z fixed to the experimental val- 
ues appropriate for the various values of x. Experimen- 
tally it is found that these parameters follow Vegard's law 
and scale linearly with x [lfj. We find that for both the 
As and P material the a and (3 frequencies vary relatively 
little with x. In Fig. [3ja) we show a weighted average of 
these calculated frequencies, F av — xFp + (1 — x)Fa s 
where Fp/F^s refers to the (average maximal and min- 
imal extremal) frequencies calculated with the atom set 
to P or As respectively. The calculated changes are much 
smaller than we observe experimentally, so it seems un- 
likely that the shrinking electron sheets can be explained 
by conventional bandstructure theory. 

An alternative is that the shrinking is driven by many- 
body effects. Recently, Ortenzi et al. [23j suggested that 
in LaFePO the shrinking of the electron and hole Fermi- 
surfaces and enhancement of effective masses can be ex- 
plained by the interaction between the electron and hole 
bands. The observed correlation of the Fermi surface ar- 
eas with the increase in effective mass and increase in T c 
would be consistent with this model; however, it remains 
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FIG. 3: color online, (a) Experimental (solid symbols) aver- 
age electron sheet frequencies (a and 0) versus P content, x. 
The data for x = are taken from Ref. 0]. The dashed lines 
show bandstructure predictions, b) The variation with x of 
the quasiparticle effective masses, m* and (c) T c after Ref. 
[l(| . The vertical dashed line marks the location of the onset 
of the appearance of magnetism at T — 0. 

to be seen if quantitative agreement can be achieved. 

In summary, we have presented dHvA data 
which shows the evolution of the Fermi surface of 
BaFe 2 (Asi_ x Pa;)2 as x is varied towards T c max . We 
find that the volume of the electron sheets (and via 
charge neutrality also the hole sheets) shrink linearly 
and the effective masses become strongly enhanced with 
decreasing x. It seems unlikely that these changes are 
a simple consequence of the one-electron bandstructure 
but instead they likely originate from many-body 
interactions. These changes may be intimately related 



to the high T c unconventional superconductivity in this 
system. 
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